Itinerancy and Hidden Order in URu2Si2 
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We argue that key characteristics of the enigmatic transition at To = 17. in URu2Si2 indi- 
cate that the hidden order is a density wave formed within a band of composite quasiparticles, 
whose detailed structure is determined by local physics. We expand on our proposal (with J. A. My- 
dosh) of the hidden order as incommensurate orbital antiferromagnetism and present experimental 
predictions to test our ideas. We then turn towards a microscopic description of orbital antiferro- 
magnetism, exploring possible particle-hole pairings within the context of a simple one-band model. 
We end with a discussion of recent high-field and thermal transport experiment, and discuss their 
implications for the nature of the hidden order. 

PACS numbers: 



I. INTRODUCTION 

The possibility of exotic particle-hole pairing leading 
to quadrupolar and orbital charge currents has been dis- 
cussed extensively in the context of the two-dimensional 
Hubbard model pi*2iiii^ More recently d-wave charge- 
density wave states, both ordered^ and fluctuatingfi have 
been proposed to explain the pseudogap phase in the 
underdoped cuprates and ground-states of doped two- 
leg Hubbard and t — J ladders. ^'^ In this paper we dis- 
cuss related anisotropic particle-hole pairing in a differ- 
ent setting, namely that of three-dimensional Fermi liq- 
uids. We believe that such pairing may occur in the 
heavy fermion metal U Ru2Si2, and here we provide the- 
oretical support for our earlier publications (with J. A. 
Mydosh) on this topic piSJiiSii^ Though the initial moti- 
vation for our orbital antiferromagnetism (OAFM) pro- 
posal in URu2Si2 was primarily experimental, here we 
observe that coexistence of large electron-electron repul- 
sion and antiferromagnetic fluctuations favours node for- 
mation in particle-hole pairing and hence the formation 
of anisotropic charge-density wave states. After present- 
ing technical details behind specific predictions for neu- 
tron scattering and for NMR, we turn towards a mi- 
croscopic description of orbital antiferromagnetism. We 
start by presenting the generalised Landau parameters 
associated with this anisotropic pairing. Next we study 
a toy model where this instability occurs. We end with 
a discussion of these results in the light of more recent 
measurements, and also suggest further experiments to 
test our ideas. 

The heavy fermion metal URu2Si2 displays a clas- 
sic second-order phase transition (see Figure 1) at 
To = n .'bK , and yet the nature of the associated or- 
der parameter remains elusive nearly two decades af- 
ter its discovery. This phase transition is characterised 
by a large entropy losatt and sharp anomalies in the 
linear^** and the nonlinear susceptibilities, ^^'^^ the ther- 
mal expansionfii and the resistivityfi^ where standard 
mean-field relations between measured thermodynamic 
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Figure 1: Schematics of the (a) specific heat anomaly. Data 
points taken from Figs. 1. of reference [13]. (b) The mea- 
sured magnetic susceptibility from reference [15] Inset show- 
ing cross-over from high temperature local Curie behaviour 
to low temperature Fermi liquid behaviour taken from Fig 2 
of reference [13]. 



quantities are satisfied^ At the transition, neutron scat- 
tering experiments observe gapped, propagating mag- 
netic excitationsSSiSiiSSiS^ that suggest the formation of a 
spin density wave. However, subsequent neutron scatter- 
ing measurements22i2^ indicate that the staggered mag- 
netic moment (toq = per U atom), is too small 
to account for the entropy loss at the transition^ which 
has been attributed to the development of an enigmatic 
hidden order. 

There is strong experimental evidence that the anti- 
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ferromagnetism and the hidden order in URu2Si2 are 
phase-separated and thus develop independently>i^ High- 
field measurenients2Si2SiSi indicate that the bulk anoma- 
lies survive up to 40 Tesla (T), while the staggered mo- 
ment is destroyedSfi by comparatively modest fields of 
15 T. Furthermore, the staggered magnetic moment 
grows linearly with pressure^^ while bulk anomalies asso- 
ciated with the hidden order remain relatively pressure- 
independent Phase separation is also indicated by 
muon spin resonance {fj,SR) experiments^i*^ The most 
direct evidence has come from recent NMR pressure- 
dependent measurements (see Figure 2)i^for T < Tq the 
existence of distinct antiferromagnetic and paramagnetic 
(hidden order) phases is clearly observed in samples with 
less than 10% of the volume magnetic {rrispin ~ O.S/is) 
at ambient pressure. The observed increase of the stag- 
gered magnetic moment with pressure^^ is then simply 
a volume-fraction effect. ^'^ The magnetic order develops 
independently from the hidden order through a first or- 
der transition^ and the associated temperature-pressure 
phase diagram has been determined using thermal expan- 



sion measurements 
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The mysterious phase transition at To has features 
that have both local and itinerant electronic natures, and 
these coexisting dual characteristics make its description 
quite challenging. For example, the development of a 
sharp propagating mode just below Tq observed by inelas- 
tic neutron scatteringS2i2ii2S emphasises the importance 
of local crystal-field excitations at the transition. Never- 
theless a purely local picture cannot provide a straight- 
forward explanation for the observed elastic anomalies^ 
near Tq that are distinct from those of typical uniaxial 
antiferromagnets'^^ both due to their (weak) magnitudes 
and due to the absence of precursor effects for T > Tq. 

The sharp mean-field nature of the phase transi- 
tion at To, together with the magnitude of the con- 
densation entropy and the observed development of 
gap in the excitation spectrum all suggest the devel- 
opment of density-wave order within a fluid of itiner- 
ant quasiparticles>ii*iSii2i2i Itineracy is implicated by the 
sharpness of the transition while gap formation and the 
large entropy of condensation speak in favour of an or- 
der parameter at a finite wavevector. However, a dis- 
senting view on this last point, involving p-wave ferro- 
magnetism, has recently been proposed;^ We note that 
within the itinerant perspective presented here, there are 
problems matching details of the excitation spectra as 
observed in inelastic neutron scattering experiments^ 
On the other hand, a purely local scenaric^^''^^ (with an- 
ticipated corrections for itinerant fermions) simply can 
not be reconciled with the almost complete quenching of 
the local moments, implicated by the paramagnetic (as 
opposed to Curie-like) susceptibihty (see inset Fig. lb.) 
and the large linear specific capacity, normally associ- 
ated with well- formed heavy electrons (Fig. la). There 
are addition inconsistencies with a local picture: for ex- 
ample, the gap A used in the local singlet schemeS^to ex- 
plain the dispersing magnetic mode has a difFerent field- 
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Figure 2: Schematic of the (a) pressure dependence of the 
ground-state staggered magnetic moment after ref. [32]. (b) 
Satellite structure in NMR after [40], taken at O.SGPa, show- 
ing the co-existence of an antiferromagnetic satellite with a 
central peak derived from the hidden order phase. 



dependence from that of the bulk A associated with ther- 
modynamic quantities.'*'^ A strict adherence to a local 
scheme requires consideration of many additional crystal- 
field levels^ evolving differently in an applied field. 

A proper theoretical description of the transition at 
To in URu2Si2 must therefore encompass both local and 
itinerant features of the problem. More specifically, the 
observed Fermi liquid properties for T > Tq (e.g. Fig. 1) 
combined with the large entropy loss and the sharp na- 
ture of the transition indicate that the underlying quasi- 
particle excitations are itinerant, presumably composite 
objects formed from the 5f spin and orbital degrees of 
freedom of the U ions. Local physics (e.g. Kondo physics, 
spin-orbit coupling, crystal-field schemes) plays a key role 
in their development. 

We have just outlined a number of general considera- 
tions that we believe are crucial features of the hidden 
order in URu2Si2- Given these criteria, we (with J. A. 
Mydosh) have proposed that it can be described by a 
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general density wave whose form factor is constrained by 
experimental observation and is ultimately determined 
by underlying local excitationsii^ We note that a num- 
ber of proposals for the hidden order that fit into this 
general framework have been niadcASiSL^SiiiiiSii^ We ar- 
gue that the large entropy loss at the transition can only 
be understood if the density-wave involves the polarisa- 
tion of a significant fraction of the quasiparticle band, a 
condition that discounts a conventional spin-density wave 
due to the small size of the observed magnetic moment. 
Taking our cue from ambient-pressure Si NMR measure- 
ments (see Figure 3) that indicate broken time-reversal 
symmetry in the hidden ordered phased we (with J. A. 
Mydosh) have proposed that UR2Si2 becomes an incom- 
mensurate orbital antiferromagnet at T = Tp with charge 
currents circulating between the uranium ions.^^ Here 
the modulation wavevector is chosen to fit the observed 
isotropic field distribution at the silicon sites. The re- 
sulting real-space fields can then be Fourier transformed 
to calculate a neutron scattering structure factor with 
a ring of possible q-vectors. Though these results have 
been presented elsewherej^''^ in this paper (Section II) we 
provide supporting technical details and further discus- 
sion. We also determine the NMR linewidths at the Ru 
sites. Detailed comparison with recent experiment puts 
constraints on the allowed incommensurate wavevectors, 
allowing us to make more specific predictions for neutron 
scattering measurements. 

In the second part of this paper, we turn towards 
an underlying microscopic picture of orbital antiferro- 
magnetism. More specifically, in Section III we explore 
particle-hole pairings in anisotropic incompressible Fermi 
liquids with specific application to URu2Si2- Next (Sec- 
tion IV) we introduce a simple t — J model with a single 
heavy band and weak antiferromagnetic spin fluctuations 
(AFMSF). We note that this particular Hamiltonian was 
originally introduced^ to describe the AFMSF-mediated 
transition in URu2Si2 at \.2K. We show that this 
same toy model also supports particle-hole pairings asso- 
ciated with incommensurate orbital antiferromagnetism 
and quadrupolar charge-density wave formation. We end 
(Section IV) the paper with a summary, and then discuss 
our results in the context of recent high-field and thermal 
transport measurements. 



II. PHENOMENOLOGY AND EXPERIMENTAL 
PREDICTIONS 

In this Section we review the experimental motivation 
for incommensurate orbital antiferromagnetism as the 
hidden order in U Ru2Si2- We develop the phenomenol- 
ogy of this proposal, independent of microscopic details. 
The magnitude and the ordering wavevector of the or- 
bital currents are fittedii to the observed isotropic field- 
distribution at the silicon sites as measured by nuclear 
magnetic resonance (NMR)iM The real-space fields pro- 
duced by the orbital charge currents at all points in the 



sample volume are then determined, and we use this in- 
formation to make specific predictions for neutron scat- 
tering structure factors and for NMR at non-silicon sites 
to test this proposal. 

A. Incommensurate Orbital Antiferromagnetism as 
the Hidden Order in URu2Si2 

We begin our phenomenological discussion by review- 
ing the case for incommensurate orbital antiferromag- 
netism as the hidden order in URu2Si2- There have 
been many proposals for the primary order parameter in 
this material)2L^2iii42iMii& and until recently it was as- 
sumed that the spin antiferromagnetism and the hidden 
order are coupled and homogenous. However pressure- 
dependent NMR measurements)^ supported by muon 
spin resonance^iiSS and thermal expansion^i data, indi- 
cate that the hidden and the magnetic orders are phase- 
separated and thus are completely independent 
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Figure 3: Temperature-dependence of the NMR line-width 
A(r), measured in Gauss, showing development of a finite 
local magnetic fields at To, after [43]. 

We believe that an important clue to the nature of the 
hidden order in URu2Si2 is provided by Si NMR mea- 
surements at ambient pressure^ that indicate that at 
T < Tq the paramagnetic (non-split) silicon NMR line- 
width develops a field- independent, isotropic component 
whose temperature-dependent magnitude is proportional 
to that of the hidden order parameter. These results 
imply an isotropic field distribution at the silicon sites 
whose root-mean square value is proportional to the hid- 
den order (ip) 

{B^{i)BP{j))^A^^H^p, (1) 

and is ~ 10 Gauss at T = 0. This field magnitude is too 
small to be explained by the observed moment22, that in- 
duces a field Bgpin = = 100 Gauss where a is the 
U — U bond length (a = 4 x 10~* cm). Furthermore 
this moment is aligned along the c-axis, and thus cannot 
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account for the isotropic nature of the local field distri- 
bution detected by NMR. These measurements indicate 
that as the hidden order develops, a static isotropic mag- 
netic field develops at each silicon site. This is strong 
evidence that the hidden order parameter breaks time- 
reversal invariance. 

The magnetic fields at the silicon nuclei have two pos- 
sible origins4i the conduction electron-spin interaction 
and the orbital shift that is due to current densities. In 
URu2Si2, the observed Knight shift indicates a strong 
Ising anisotropy of the conduction electron fluid along the 
c axis; therefore the electron-spin coupling is unlikely to 
be responsible for the measured isotropic field distribu- 
tion at the Si sites. It thus seems natural that these local 
fields are produced by orbital currents that develop at Tq, 
and thus we attribute the observed isotropic linewidth to 
the orbital shift. 

It is this line of reasoning that led us (with J. A. My- 
dosh) to propose^-'^ that URu2Si2 is an incommensurate 
orbital antiferromagnet at T = To with charge currents 
circulating between the uranium ions. The planar tetrag- 
onal structure of URu2Si2 presents a natural setting for 
an anisotropic charge instability of this type. We can 
estimate the local fields at the silicon sites that are pro- 
duced by the orbital currents. On dimensional grounds, 
the current along the uranium-uranium bond is given by 
/ ^ where A is the gap associated with the formation 
of the hidden order at Tq] we note that this expression 
also emerges from an analysis of the Hubbard model i 
If this orbital charge current is flowing around a ura- 
nium plaquette of side length a, then the magnetic fleld 
produced at a height a above it is given by Ampere's 
Law to be i? « = IIG, in good agreement with 

the observed field strength;— here we have used the ex- 
perimental valuai^ A — 110 A'. Note that the resulting 
orbital moment, moAFM — 0.02/iB {moAFM — lo^), is 
comparable to the effective spin moment at ambient pres- 
sure. We emphasise that an orbital moment produces a 
field an order of magnitude less than that associated with 
a spin moment of the same value; the low field strengths 
observed at the silicon sites are quantitatively consistent 
with our proposal that they originate from charge cur- 
rents. 

This orbital moment, moAFM — 0.02/1^, can also ac- 
count for the entropy loss at the transition. We empha- 
sise that its large value suggests that the amplitude of 
any proposed density-wave must be a significant fraction 
of its maximally allowed value, and will proceed to show 
that this is the case for the OAFM. In a metal the change 
in the entropy is given by AS" — A7„To where A7„ is 
the change in the linear specific heat coefficient resulting 
from the gapping of the Fermi surface. A7„ is inversely 
proportional to the Fermi energy e^? of the gapped Fermi 
surface, so that in general the change in entropy per unit 
cell is given by A5 = ^ (^f^)- Since the transi- 
tion at To is mean-field in nature^ we have A ~ To so 
that A5 ~ — . Now we recall that the orbital magnetic 



moment is 



moAFM = la" = (J) a^A « 0.02/is (2) 
such that it is saturated when A ~ ej? 



'■OAFM 



analogous to the saturation value of the electron spin 
/is — (f ) iQ^i/ where ag and e[j are the Bohr radius and 
the energy of the Hydrogen atom respectively; here we 



have used 



10^ and ^ 



M* 



lO"'^ where Mh 



"0 <^H J" 

and M* refer to the mass of hydrogen and of URu2Si2 re- 
spectively. Then the change in entropy at the transition 
(AiS ^ due to the development of orbital antiferro- 
magnetism can be expressed as 



A5o 



AFM 



( moAFM 



OAFM 



0.02 



"-OAFM 



0.2 (4) 



which is a number (0.2 — 0.3 In 2) in good agreement 
with experiment pi^ We also note that the critical field 
for suppressing the thermodynamic anomalies is distinct 

from its spin counterpart: the ratio 



10 is qualitatively consistent with the observed critical 
field associated with the destruction of hidden order 
We emphasise that the sizable entropy loss associated 
with the development of orbital antiferromagnetism in 
URu2Si2 is a direct consequence of its renormalised elec- 
tron mass oc More generally the orbital mo- 
ment is a larger fraction of its saturation value than is 
its spin counterpart, and this leads to the large entropy 
loss. 

Orbital antiferromagnetism can therefore account for 
the local field magnitudes at the silicon ions and for the 
large entropy loss at the transition. Our next step is 
to tune the ordering wavevector to fit the isotropic dis- 
tribution at these sites and then to determine the real- 
space fields throughout the sample volume. This can then 
be Fourier transformed to make predictions for neutron 
scattering. We note that it has been suggested'^^ that 
the isotropic nature of the field distributions at the sil- 
icon sites may be due to impurity-broadening. Though 
disorder is certainly present in these samples, we believe 
that the incommensurate nature of the density wave is 
the origin of this isotropy. Towards proving this point, 
we have determined the anisotropic field distributions at 
non-silicon sites; their observation via NMR would cer- 
tainly not be possible if there were significant disorder- 
smearing. 

Before proceeding with this program, let us com- 
ment brieffy on the current experimental situation re- 
garding the proposal of incommensurate orbital anti- 
ferromagnetism in URu2Si2- We admit that our pro- 
posal is closely linked to the ambient-pressure NMR 
experimentSjM which are the only direct evidence of bro- 
ken time-reversal symmetry in the hidden ordered phase 
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and have not been reproduced by other groups. We note 
that muon spin resonance measurement o'^^''^^ support the 
emergence of local fields with the same temperature- 
dependence as that associated with NMR, but their over- 
all amplitudes are two orders of magnitude less than 
that seen in the NMR measurements. This is a point 
to which we return in the discussion. Although incom- 
mensurate peaks have been seen in inelastic neutron scat- 
tering measurements these are due to excitations 
above the partly gapped Fermi surface and are not di- 
rectly related to the orbital antiferromagnetism. Cur- 
rent experimental resolution for elastic scattering - a di- 
rect probe of the incommensurate orbital antiferromag- 
netic order - is not yet good enough to confirm or deny 
the OAFM scenario. Here we present technical support 
for previous predictions for neutron structure factors/^ 
while also making specific suggestions for measurements 
where the signal should be sufficiently strong to be ob- 
served practically. 



as shown in Fig. 01(a). The circulating current around 
plaquette j is then taken to be 



/c(Xj) = /oe^Q-^J +H.C 



(6) 



Using Ampere's law, link 1-2 will will produce a contri- 
bution to the vector potential given by 

A-(x) = iV r^''d.'MM^ii, (7) 



where Xi2 is the unit vector pointing along the bond from 
1 to 2. Writing x'j as 

x', =x5.^'+w(x(2)-x«), 
where < w < 1 defines the position along the link, we 
have 



B. Predictions for Neutron Scattering 

In order to calculate the neutron cross section for scat- 
tering by incommensurate orbital antiferromagnetic or- 
der, we use the Born scattering formula, 



dw 



/c(Xj-) Xi2 



|x - {xf ^ + i(;(x(2) -x(i))}|' 



(8) 

for the vector potential where a is the U-U bond length 
in the ab plane. 



da 
dh 



9Ne 



|B(q)p = r^^(q), 



(5) 



where is the neutron gyromagnetic ratio, q the scat- 
tering wavevector of the neutrons, |B(q)|^ is the struc- 
ture factor of the magnetic fields produced by the orbital 
currents and 5(q) = |B(q)p/(47r^B)^ is the structure 
factor measured in units of the Bohr magneton {(J-b)- 

We shall compute the magnetic field as the curl of the 
vector potential, B(x) = V x A. The procedure will be to 
compute the vector potential produced by the circulating 
current around a given plaquette. We shall denote the co- 
ordinate of the centre of plaquette j by Xj . The corners 

of this plaquette are located at sites x^'^'', (r = 1, 4) where 



X . - 



(r = l,4). 




(b) 



T iQX 



i n e ^ ^ 




Figure 4: (a) Labelling of sites around a single plaquette. 
(b) Schematic of circulating currents /c(Xj) = /oexp[iQ • 
Xj] flowing in the uranium plaquettes in the ab plane. The 
plaquettes are labelled by the coordinates of their centre Xj. 



We now compute B^-^ = V x A^^, and take the Fourier transform to obtain 

B-M = ^T./y-'l"'^ 'c(^, ^ I. - (X, + X... ! Mx<" - X..) 



)}l 



c 



E/c(Xj) xi2 X q / dw d 



|x - {xj + x(i) + ii;(x(2) -x(i))}|' 



(9) 



Using 

/I 47r 
d3g-.qx.^ = — e-^i"", 

|x - a| 



6 



we obtain 



B^^(q) 



iAna 



V/c(Xj)ii2 xq f dw exp[-iq.(xj +x(i) +«;(x(2) -x(i)))] 

"'0 



Ana 



q.(x(2) -x(i)) 

J 



(10) 



where we have replaced /c(Xj) = /oe**^ -^^ and 



Xl2 



(1) 



F^^(q) = i^(e-»q"'^ _ ) 



q ■ xi2 



(11) 



is the form- factor associated with hnk 1 — 2 in the plaquette centred about Xj. To sum over ah of the hnks around 
the plaquette, we must add together the form factors 



F(q) = Fi2(q)+F2^(q)+F-^4(q)+F^i(q 

X 

q.x 

^_g-iq-(-x+y)a/2 _ g-iq.(x+y)a/2| _ x ^ y 



.|g«q-(*+y)a/2 _ g-iq.(x-y)a/2 



4 sin 



2 



qytt 
2 



_y ^ 

qy qx 



(12) 



We let Q be the wavevector for the orbital order so that /(xj) = Iq exp[— zQ.Xj]. 

Replacing F-'^^ — > F in Eq. 1101 we we obtain the complete Fourier transform of the magnetic field: 



B(q) = exp[«(Q - q)-Xi] sin sin 



qy_a 
2 



q • y qx 



X q. 



(13) 



The U sites Xj can be written as 



X, = a (ji, j2,0) + I (0,0, J3) + i(l - (-1)^^) (|, ^,0), 

(14) 

where c is the separation between even or odd num- 
bered U planes. The unit cell has lattice vectors 
(a, 0, 0), (0, a, 0), (0, 0, c). For an isotropic distribution of 
magnetic fields at the Si sites, we can reasonably expect 
Q to be staggered between successive U layers. We per- 
mit Q to be incommensurate in the a — b plane: 



Q = (Q.,Q„0) 



27r 



(0,0,1) 



(15) 



Summing over the lattice sites in Ea. (|13(l we find 
Svr/o 



B 



OAFM 



(q) 



q'^c 



q,Q+G 



1 



JG-{a/2,a/2,c/2) 



X sm 



m(— j sn 



qya 



y_ 

q.x q.y 



where G = 27r [ni /a, 712/0, ns/c] is a reciprocal lattice 
vector. 



Equation (|16l) should be contrasted with the corre- 
sponding expression if the order parameter were a spin 
density wave instead of an orbital antiferromagnet: 



B 



SDW 



(q) 



An 



q.Q+G 



1 +e»G-(a/2,a/2,c/2) 



{q X (M X q)} . 



(17) 



We note that a major difference between the two cases is 
that BoAFM(q) decreases rapidly as q~'^ while Bs_Diy(q) 
is constant. This makes OAFM much harder to detect 
in neutron scattering experiments than its SDW coun- 
terpart. Second, the term (q x (M x q)) in Boa_fm(?) 
indicates that scattering is suppressed for q = Q since 
for an SDW along the c axis, M || Q = ^(0, 0, 1). here 
is no such term in ^oafm{q)- Thus the presence of a fi- 
nite scattering amplitude at this particular wavevector in 
URu2Si2 would be a "smoking gun" confirmation of in- 
commensurate orbital antiferromagnetism as the hidden 
order. 

Next we turn to obtaining the structure factor |B(q)|2. 
Neutrons couple to the orbital currents via their magnetic 
moment (/ijy = SatM-bS) as E = — /x^.B. For incoherent 
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neutrons, |B(q)p is the modulus squared of Q averaged over the orientation, the neutrons. Thus 



^(q) 



|B(q)|2 fNW 



Sq.Q+G [jo 



1 + cos[7r(ni +n2+ ris)] 



r 



'Qxa- 




Qya- 


I 2 \ 


jo 


I 2 \ 

















(18) 



where joix) = and N is the number of U sites. From 
this expression, we find that the maximum scattering in- 
tensity is predictedii to he in a ring Q = Qq + q oi 
radius \q\ ~ 0.2 centred on the wavevector Qo = (001) 
where q'hes in the a — b plane. Once again, we emphasise 
that scattering in the vicinity of Qq is forbidden for the 
case of ordered spins along the c-axis; thus the observed 
presence of neutron scattering intensity at this particu- 
lar wavevector would be a "smoking gun" confirmation 
of orbital antiferromagnetism as the hidden order. 

In general the structure factor can be written as a prod- 
uct 



S{q) = f{q)9{q) 



(19) 



where g{q) is a function periodic in the reciprocal lattice 
vector but /(g) is not. For the case of orbital antifer- 
romagnetism, the calculated structure factor yields an 



J 



asymptotic form for the form factor f{q) ^ This 
power-law decay of the intensity peaks is due to the ex- 
tended nature of the scattering source in contrast to the 
exponentially decaying structures observed for point-like 
spin antiferromagnetism. 



It is tempting to state that such power-law peaks will 
be a clear signature of orbiting charge currents, but we 
still need to determine whether the overall intensities are 
observable. We can estimate the strength of the pre- 
dicted OAFM neutron signal compared to that associated 
with spin magnetism at ambient pressure. Our calcula- 
tions indicate that a fifth of the total integrated weight of 
S{q) (TIWSQ) resides in the first Brillouin zone for the 
OAFM. Using the sum rule that relates the total ISWQ 
(integrated weight of S{q) to the square of the moment, 
we have 



{IWSQ)bzi 



\{TIWSQ)oAFM = \{moAFMf' 



1 

500 



1 
500 



{TIWSQ)s 



(20) 
(21) 



r 



where we have used tuoafm = 0-2^b and rrispin = 
0.3/iB. Since the magnetic region occupies roughly a 
tenth of the sample at ambient pressure we then write 



{IWSQ)bzi^-^ MeasnvediT IWSQ)sp^n (22) 

which indicates that the scattering peaks in the first Bril- 
louin zone due to orbital ordering should have roughly 
1/50 the intensity of the analogous spin peaks at am- 
bient pressure. There have been two exploratory neu- 
tron studies^2i^ but neither was conclusive due to is- 
sues of resolution. In particular the more recent elas- 
tic measurements^ were not performed at the predicted 
wavevector Qp = (rp, Tp, 1) where there should be no 
dipole scattering; please recall that here the form factor 



~ ((f X m) and m is aligned with the c-axis. More specifi- 
cally the scattering intensity should be a factor of twenty 
higher than at Qe = (1 + Tx, Ty, 0) where the experiments 
were performed, and the experimental resolution should 
be good enough then to prove/refute the orbital antifer- 
romagnetism proposal. 



Nuclear Magnetic Resonance Linewidth at the 
Si and the Ru Sites 



Nuclear Magnetic Resonance (NMR) is a local probe 
of the strength and the local distribution of the magnetic 
field distribution in the material. We use experimental 
NMR results to determine the ordering wavevector as- 
sociated with the orbital antiferromagnetism, which can 
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then be included in the structure factor calculated above. 
Thus neutron scattering and NMR are complementary. 
Eq.(jHl gives the vector potential at a point x due to a 
current in link (12) of a plaquette centred at Xj. Contri- 
butions from other links in the plaquette may be similarly 
written out (see Fig0J. The magnetic field at any point 
X can be obtained using B = V x A, where A is the total 
vector potential obtained by summing contributions from 
all links and plaquettes. We give detailed expressions for 
A in the appendix. 

For the sake of completeness, we reviewiS^ii our ar- 
guments regarding the Si NMR measurements''^ and the 
ordering wavevector of the orbital antiferromagnetism. 
We note that the silicon atoms in URu2Si2 are located 
at low-symmetry sites above and below the uranium pla- 
quettes, so that the fields there do not cancel. Therefore 
the proposed OAFM must have an incommensurate Q 
(tt, tt) in order to produce isotropic field distributions at 
the silicon sites. If the order parameter in the hidden 
order phase is OAFM, then such a magnetic field distri- 
bution at the Si sites would be possible if the wavevector 
for orbital ordering were incommensurate^Siii 

Q = —(0.22 COS0, 0.22 sin 0,a/c). (23) 
a 

Fig. (0 shows the distribution of the magnetic field lines 
about the ab plane for an incommensurate Q correspond- 
ing to (/) = 7r/4 in Ea. lpS|l . and viewed in the [DID] direc- 
tion. 



0.25 0.5 




0.25 0.5 

Qx 

1 



Figure 6: Anisotropy of magnetic field distribution 
at the Si sites. The anisotropy is found to vanish 
on a ring of wavevectors approximately given by Q = 
^(0.22 cos 0, 0.22 sin 0, a/c) . 



t ^^^^^^^^^^^^^^^^^^ ^ ^ ^ 

Figure 5: The distribution of the magnetic field lines about 
the ab plane for an incommensurate Q = 2lL (0.16, 0.16, a/c), 
and viewed in the [010] direction. The black circles represent 
Uranium atoms in the ab plane. 

A convenient definition of the anisotropy in the mag- 
netic field at a given site is 

C= |(i?±-S||)/(i?±+B||)|, 

Fig.® shows the anisotropy as a function of the Q vec- 
tor. While the field distribution at the Si sites is isotropic, 
that need not be the case at other sites such as Ru; fur- 
thermore the anisotropic nature of the field distribution 
at the Ru sites would indicate that disorder-averaging is 
not at play here. If we take as the origin any uranium 
atom in the lattice, we find the Ru sites at coordinates 

Xru = |(*-J + l>*+J-0) + ^(0-0,/c+l/2), (24) 

where i,j,k are integers. FigEl shows the anisotropy of 
the magnetic field distribution at the Ru sites. 



Recent Ru NMR measurements^ report a local mag- 
netic field anisotropy of around 0.3. Values of Q deduced 
from our OAFM model using the Ru NMR data should 
of course be consistent with Si NMR. The anisotropy of 
the magnetic field at the Ru sites calculated from our 
model shows strong variations as the orientation of the 
incommensurate wavevector given in Ea. H23() is varied. 
Anisotropy of field at the Ru sites for OAFM ordering 
wavevectors given by Ea. (|23|l varies from about 0.7 along 
the (j> = 0, 7r/2 directions to nearly unity along (j) = n/A. 
Thus the most likely incommensurate wavevector Q for 
OAFM lies close to the (j) — 0,7t, ±7r/2 directions. 

Neutron scattering measurements^ show enhanced 
scattering for T > Tq at the incommensurate wavevec- 
tors 

Qexp ~ (27r/a)(ni -|- 0.4cos(/), n2 + 0.4cos(/), 713), (25) 

where ni + n2 + is an odd integer. Below Tq, the 
ring of excitations seems to collapse toward the x and 
y directions, decreasing in intensity. The structure- 
factor predicted in Ea. (|18|) could not be verified/refuted 
due to issues of resolution."*^ According to Ea. (|18|l . the 
structure factor measured near Qexp = (27r/a)(1.4, 0, 0), 
as was done in the most recent experiments^ has a 
scattering intensity that is smaller than that at Q = 
(27r/a)(0.4, 0, a/c) by a factor of more than five. In 
an earlier experiment enhanced scattering was ob- 
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Figure 7: The anisotropy of magnetic field distribution 
at the Ru sites for orbital wavevector of the form Q — 
(Qi, Qa, Stt/c). Darker shades indicate lower anisotropy. 

served at Q = (27r/a)(1.4, 0, 0) above the transition tem- 
perature Tq. The scattering intensity was sharply en- 
hanced for T < To, and furthermore, the scattering 
linewidth decreased to resolution-limited values. More 
work is needed to verify whether the incommensurate 
peak observed in neutron scattering measurements is re- 
lated to Eq.|j23|l deduced from Si and Ru NMR data 
using our model of orbital antiferromagnetism, and we 
strongly suggest elastic neutron scattering measurements 
at the wavevector predicted to have the greatest intensity 
(Q = (27r/a)(0.4, 0, a/c)) to test OAFM as hidden order. 

III. TOWARDS A MICROSCOPIC 
DESCRIPTION OF THE HIDDEN ORDER 

We now turn to a more microscopic approach to the 
hidden order. As we have already noted, a proper the- 
oretical description of URu2Si2 must encompass both 
local and itinerant features of the problem. A gen- 
eral duality scheme for heavy electron systems has been 
proposed^ In this model, the itinerant excitations are 
constructed from the low-lying crystal-field multiplets of 
the uranium atom. The quasiparticles associated with 
the heavy Fermi liquid are composite objects formed from 
the localised orbital and spin degrees of freedom of the U 
ions and the conduction electron fields. The phase tran- 
sition in this model is then a Fermi-surface instability of 
these composite itinerant f-electrons. This approach has 



been adapted^ to describe the coexistence of hidden or- 
der with a small moment in URu2Si2- With the more 
recent understanding that the hidden ordered phase does 
not contain a staggered magnetisation, we have revisited 
this duality scheme and, guided by experiment, now 
discuss its implications for the nature of the mysterious 
order that develops at T ~ Tq. 



A. Possible symmetries for particle-hole pairing 

We begin with the assumption that all the excitations 
of URu2Si2 that condense into the hidden ordered state 
are of itinerant character. More specifically, we will as- 
sume that all of the system's local physics (e.g. local 
moment character of the f-electrons) has been absorbed 
into the formation of composite quasiparticles. Given 
this premise, it then follows that key aspects of the (hid- 
den) order parameter will be expressed through its ma- 
trix elements between quasiparticle states. If we denote 
it by the operator , then its general matrix element 
between quasiparticle states is 

(k + Q/2, a|^P|k - Q/2, a') = A^' (Q) (26) 

where Q is the ordering wave- vector and jkcr) is the quasi- 
particle state of momentum k. Microscopically we would 
have to characterise ^ in terms of the detailed crystal- 
field split states of the U ion, but for the purposes of 
characterising the phase transition, quasiparticle matrix 
elements should suffice. Within the Hilbert space of the 
mobile f-electrons, the order parameter can then be writ- 
ten 

#^Ar'(Q)c^Q/2^,Ck.Q/2,.'. (27) 

where (Q) is a general function of spin and momen- 
tum. 

We are therefore considering a class of density waves 
with the most general pairing in the particle-hole chan- 
nel characterised by (Q). We now categorise the 
possible particle-hole pairingsiS in URu2Si2- Assum- 
ing that the hidden order develops between the uranium 
atoms in each basal plane, we restrict our attention to 
nearest-neighbour pairings on a two-dimensional square 
lattice, and display the five resulting possibilities in Ta- 
ble in Ea. H27|l . We emphasise that each of these pair- 
ing choices will partially gap the Fermi surface, account- 
ing for the large entropy loss and the observed anoma- 
lies in several bulk quantities .^^ In conventional charge- 
and spin-density waves (CDWs and SDWs respectively), 
the quantity ylk(Q) is an isotropic function of momen- 
tum. However in more general cases Ak(Q) will develop a 
nodal structure which leads to anisotropy (Table P) that 
is favoured by strong Coulomb interaction, as we shall 
discuss in the next section. 
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Name 


^r'(Q) 


T-invariance 


Local Fields 


SDW (isotropic spin density wave) 


(7 


no 


yes 


CDW (isotropic charge density wave) 


const. 


yes 


no 


d-SDW 


(7 {cos{kxa) — cos(fcya)) 


no 


no 


q-CDW (quadrupolar) 


cos{kxa) — cos{kya) 


yes 


no 


OAFM (orbital antiferromagnet) 


i (sin(A;a;a) — sin(fcj;a)) 


no 


yes 



Table I: Possible symmetries for particle-hole pairing 



B. General Discussion of Anisotropic Charge 
Instabilities in Fermi Liquids 

At low temperatures, heavy electron materials form 
almost incompressible Landau Fermi liquids in which 
the residual interactions between heavy quasiparticles are 
driven by strong, low-lying antiferromagnetic spin fluctu- 
ations. This harshly renormalised electronic environment 
is conducive to the development of instabilities in which 
electrons or holes form bound-states that contain nodes 
in their pair wavefunction. 

Such arguments are well established in the context of 
anisotropic Cooper pairing^^S Here we extend these 
ideas, arguing that an almost incompressible Fermi liq- 
uid is highly susceptible to the formation of anisotropic 
density waves, where the staggered electron- hole conden- 
sate contains a node in the pair wavefunction. This 
issue first arose in the context of orbital ordering in 
cuprate superconductors^. Here it has been empha- 
sised that strong Coulomb interactions suppress electron- 
hole bound-state formation in CDWs, unless the bound- 
state contains a nodei& Heavy electron fluids provide a 
unique opportunity to apply these arguments to three- 
dimensional systems. Furthermore there is no contro- 
versy associated with the Landau-Fermi liquid of their 
normal states, a situation in distinct contrast to the sit- 
uation in the cuprates. 

In a heavy electron fluid, the density of states is 
severely renormalised so that the ratio of the quasiparti- 
cle and the bare band-structure density of states 



A^*(0) 



1 



is typically at least a factor of ten. In these systems 
the magnetic susceptibility, given in Landau Fermi liquid 
theory by 



l + Fi' 



is weakly enhanced. By contrast, the charge suscepti- 
bility is severely depressed by strong coulomb interac- 
tions and is essentially given by the unrenormalised band- 
structure value 



jY*(o) 



which is why the fluid is characterised as "almost incom- 
pressible". It is this basic effect that rules out the for- 
mation of isotropic charge density wave order and s-wave 
superconductivity. 

Response functions that contain an anisotropic form 
factor are unaffected by the strong Coulomb interactions. 
The key point here is that the strong interaction effects 
are local and thus they do not affect the higher Lan- 
dau parameters, due to the nodes in the corresponding 
spherical harmonics. For example if we consider a "chem- 
ical potential" which couples anisotropically to the Fermi 
surface in the 1-th angular momentum channel, then the 
corresponding susceptibility is given by 



X. 



(0 



N*{{)) 



1 



N*{Q) 



N{Q) 



provided the higher Landau parameters are not much 
larger than unity. From this discussion, we see that large 
mass renormalisation and strong Coulomb repulsion sup- 
presses isotropic CDW formation but that analogueous 
instabilities can form in higher angular momentum chan- 
nels. 



C. Anisotropic Pairings: The Contenders 

We have just argued that the large Coulomb repulsion 
between the heavy fermion quasiparticles (incompress- 
ibility) in URu2Si2 discourages isotropic pairing in the 
CDW channel. This expectation is confirmed by exper- 
iment, for charge density wave formation is expected to 
produce a lattice distortion, yet none is observed to de- 
velop URu2Si2 below the 17K phase transition. Simi- 
larly neutron scattering is inconsistent with the presence 
of an isotropic spin density wave in the hidden ordered 
phaseiSSiSi Thus, mainly due to the incompressibility of 
the heavy Fermi liquid, we are left with three remaining 
anisotropic particle-hole pairing states (see Table QJ. 

The possibility of d-spin density waves as the hidden 
order in URu2Si2 has been raised by several authorsiiSiS 
In a Stoner analysis, d-SDWs require ferromagnetic ex- 
change interactions of neighbouring spins. In particu- 
lar, for antiferromagnetic interactions, a Stoner analy- 
sis reveals that the d-SDW has a lower transition tem- 
perature than competing quadrupolar CDW (q-CDW) 
or spin-density waves.^^ Thus a d-SDW scenario favours 
ferromagnetic fluctuations in U Ru2Si2] by contrast, its 
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transition at T* = 1.2K to a d-wave superconductor in- 
dicates the importance of antiferromagnetic fluctuations 
at T > T*. 

Before discussing the two remaining options presented 
within the framework of Table ^ we want to mention 
two recent proposals for the hidden order parameter 
that both lead to quasiparticle matrix elements similar 
to those of a higher-order SDW. In the first one, the 
authors^ argue that consistency with experiment can be 
maintained for an SDW that develops predominantly in 
the p- or s- bands whose neutron form-factor at the Bragg 
peaks is significantly smaller than that of f-electrons. 
Here the key conceptual difficulty is that the matrix ele- 
ment of the order parameter in the f-bands would have to 
be small; yet the large entropy of condensation observed 
at T = To is almost certainly associated with these same 
f-electrons. It has also been suggested^ the hidden order 
results from octupolar crystal-field states. In the quasi- 
particle basis, such an order parameter behaves like a 
spin-density wave with a small g— factor. At present, 
the viability of this approach awaits more detailed pre- 
dictions regarding the magnetic distributions within the 
sample that then, like for the OAFM scenario, could be 
tested by NMR and neutron measurements. 




/(k) = cos(A:^) - cos(fcj^) /(k) = i[sm(fc^) - sin(ky)] 



Figure 8: (a) Incommensurate quadrupolar density wave 
(qCDW). In two dimensions, the form factor cos(fca;) — cos(fcj;) 
leads to a an incommensurate density wave with a quadrupo- 
lar charge distribution, the CDW analogue of a d-wave su- 
perconductor, (b) Incommensurate orbital antiferromagnet. 
Here currents circulate around square plaquettes defined by 
nearest-neighbour uranium ions. 

Returning to the table of possible pairing symmetries 
(Table P) we therefore have two remaining options: the 
quadrupolar charge density wava^ (Fig. |Hl(a)) and the 
orbital antiferromagnet (Fig. |Sl(b))>i2iii, where both sce- 
narios are consistent with our picture of URu2Si2 as 
an incompressible Fermi liquid with strong antiferromag- 
netic fluctuations. Each order parameter has nodes, so 
that neither couples directly to the local charge density. 
Furthermore both incommensurate density waves couples 
weakly to uniform strain, and thus are both consistent 
with the observed insensitivity^ of the elastic response 
at Tq. Recent uniaxial stress measurements suggests that 
the hidden order is sensitive to the presence of local 
tetragonal symmetry^S a feature that can be explained 
within both frameworks for completely different reasons. 



In the orbital antiferromagnet the currents are equal in 
each basal direction^ ^, whereas within the quadrupolar 
scenario it is known that some of the crystal-field states 
with tetragonal symmetry are quadrupolar4S Unfortu- 
nately the diamagnetic response cannot be used to dis- 
criminate between these two scenario, as the contribu- 
tion from orbital antiferromagnetism is small compared 
to that associated with the gapping of the Fermi surface 

At present, the key factor distinguishing the orbital an- 
tiferromagnet from the quadrupolar charge-density wave 
scenarios is the absence or present of time- reversal break- 
ing. Because the local field distributions and strengths 
measured by NMR have not yet been observed by other 
methods, there is still uncertainty about these results. 
We note that it has been arguedii that the observation 
of a stress-induced moment^ implies that the hidden- 
order breaks time-reversal symmetry; much as we would 
like to believe this, we note that this result can be at- 
tributed to a volume-fraction effect and thus is incon- 
clusive. Both the quadrupolar charge density wave and 
the orbital antiferromagnet have nodes in their respective 
gaps, which should in principle be observable via photoe- 
mission and/or scanning tunnelling microscopy, though 
issues associated with the nature of the surface of this 
material remain to be resolved. However the quadrupo- 
lar charge density wave is not expected to lead to mag- 
netic neutron scattering, and therefore detailed elastic 
measurements are critical for resolving the nature of the 
hidden order parameter. 

IV. TOY MODEL FOR ANISOTROPIC 
PARTICLE-HOLE PAIRING 

Next we explore a simple t — J model for heavy elec- 
trons with antiferromagnetic spin fluctuations, and ex- 
plore different orderings. We are motivated by exper- 
iment in our choice of the model. URu2Si2 undergoes 
a phase transition to a d-wave superconducting state at 
To = 0.8K, and the pairing is understood to be medi- 
ated by antiferromagnetic spin fluctuations. The same 
t — J model also encompasses orbital antiferromagnetism, 
quadrupolar CDW, and isotropic SDW. 

We consider a simplified model for the heavy Fermi 
liquid, described hy H = Hq + Hj, where 

k 

describes the band of heavy electrons and 

iJz=^J(q)S(q).S(-q), (28) 
q 

is the interaction between them. Here, S(q) = 
k^ti+qa'^ afjCkp IS the Fouricr transform of the local spin 
operator. . In this simplified model, we consider the in- 
dices (7 to represent the pseudo-spin indices of the spin- 
orbit-coupled, heavy-electron band. We recall that we 
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are working in an itinerant basis where the local physics 
(e.g. spin-orbit coupling) is absorbed into the compos- 
ite quasiparticle states. Using the completeness relation 
(Tap ■ S^s + SapS^jj — 2dariS-yp wc may rewrite this inter- 
action as 



H, 



ij, a a' 



1 



Here we have rewritten the electron operators in a local 
basis, so that Cja = c\^^e^^'^^ is the electron cre- 

ation operator at site j, N is the number of sites in the 
lattice and Jy = jf'^qJ{(i)e^'^'^^^~^^ is the spin inter- 
action between sites i and j. We shall ignore the second 
term, which involves the heavily suppressed fluctuations 
in quasiparticle occupation at each site. The first term 
can be decoupled as 

= E^(l^- P)<k+c.,k_4^p_c,',p+,(29) 



q.k.p.crcr' 



where k± = k ± iq, etc. The interaction potential J(q) 
can be expanded into partial waves, 



Vi=2 dxxPi{l- 2x^) J{2px) 



(30) 



where x — sin(0/2) and p « pp- Wc require that the 
I = (isotropic) component be large and negative, re- 
flecting strong on-site quasiparticle repulsion. This has 
the effect of suppressing isotropic particle-hole pairing. 
However this potential VJ, for / > 0, could be attractive, 
which would then favour particle-hole pairing in higher 
angular momentum channels. Such higher angular mo- 
mentum components are present due to anisotropy of the 
interaction J(q), which occurs at sufficiently large values 
of q where the underlying symmetry of the crystal be- 
comes important. 

For the purposes of a toy model, we shall assume in 
URu2Si2, nearest neighbour antiferromagnetic spin fluc- 
tuations (AFMSF) predominate, so that 



J(q) «2Ji7i 



(31) 



where the form factor 7q = cos(ga;a)-t-cos((7j,a). With this 
approximation, the interaction in Ea. (|29|l is separable: 



qeiBZ,fc,p;r=l,4 



(7;^Pp(q))H^Pk(q), (32) 



are form factors that transform under the point-group 
symmetry of the lattice. Since /9k(q) = Pk(— q), the 
quantity inside the summation is symmetric under q — > 
— q, and so, by doubling the prefactor and restricting the 
sum over q to one-half the Brillouin zone, we assure that 
every term in the q sum is independent. 

This interaction is attractive and of equal magnitude 
in the four anisotropic channels. 7^, 7^, 7^'^ have s- 
like, d-like and p-like symmetry respectively. Notice that 
bond-variables ^„{c\„Cj<j) are invariant under time re- 



versal, Y.a{4aCja) = Y.a{'^ja<^^'y)* ^.nd the imaginary 
pre-factors in 7^''* have been chosen so that the form- 
factors respect this symmetry, i.e 7^ = (7?^^.)*. 

By carrying out a "Hubbard Stratonovich" decoupling 
of i//, wc obtain 



Hi 



E 



A^7^Pk(q) + A^(7D>lc(q) 



N 
2Ji 



BZ,k:r=1.4 



E 



(35) 



qeiBZ; r=i,4 



Now the mean-field solution to this expression is deter- 
mined by the saddle-point condition 



^q--#E(7£)*(^'^(-q)) 



N 



(36) 



In general, the density wave will condense at a primary 
wavevector q = Q. For a realistic model, Q may well be 
incommensurate, in which case, it will be accompanied 
by a family of corresponding Q' that form a "star" of 
q-vectors under the point group. There will in general 
also be higher harmonics of Q. To illustrate the key 
ideas however, we shall assume a simple model in which 
a single Q dominates the density wave, i.e. 



A^ = A^ 
q 



A^<5, 



q -Q 



For this discussion, we shall also assume that the Fermi 
surface is "almost nested" , so that the Fermi surface can 
be divided into two equal parts or reduced Brillouin zones 
(RBZ): region I in which Ick-iql ^ kk-i-iql and region 
II in which Ifk-iql ^ kk+iql- perfectly nested Fermi 
surfaces Ck = — Ck+q are perfectly degenerate. For a 
square lattice and Q = (tt, tt) the reduced Brillouin zone 
is the diamond-shaped region bounded by — tt < ky < 



TT, — TT - 



\ky \ < fcj, < TT 



The mean-field Hamiltonian is then 



where 



Pk(q) 



E '^k+iqtr'^k-iqfT 



are the particle-hole operators and 



(33) 



H 



MFT 



E ^'-k+'^k- 

r=i,4 



Ji 



ek+ 


Ak' 


(::) 


Ak 


fk- 





(37) 



1.2 
7k' 

3.4 
7k 



cos(fca;a) ± cos{kya) 
i{sui{kxa) ± sin(A:j,a)) 



(34) 



where k=^ = k ± Q/2 and Ak = A^7k- Here, Q is 
the wave vector for particle-hole pairing. Diagonalising 
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the electronic part of the total Hamiltonian yields two 
bands, 



(±) 



1 



±\/4(ek+iQ-ek-iQ)' + |Ak| 



(38) 



The mean field solution for pairing density (pp(Q)) in 
Ea. (|S5|l is obtained by setting the variation of the free 
energy F, 



RBZ 



F ^ -TY, ln[(l + exp(-/34+)))(l + eM~PE^^^^))] 



Ji 



r=i,4 



(39) 



with respect to A'" to zero. This yields the gap equation, 

RBZ ^'rj^i-)-^ _ f(T^i+)^ 

(40) 



AT 



1 



N ^ 
k.r 



^ A,(,r). /K-')-/(^r') 



4*' - eV 



In the special case where condensation occurs in a single 
channel F = Fq, this simplifies to 



Ji 



RBZ 

- = -y 



7k 



ro,2/(4'V/(4+^) 



4+) -4-) 



(41) 



At T = To, Equation H4U|) is essentially a Stoner criterion 
'^iXoV'(O) = 1 where 



RBZ 



xo^(q) 



N 



0|2 /(gk--q/2) - /(£k++q/2) 
£k++q/2 — ek-_q/2 



(42) 



is the susceptibility associated with the hidden order pa- 
rameter ^/;, measured at a wave vector q + Q. 

Without details of the band-structure wc can not pre- 
dict which of the four order parameters will dominate. 
Some general comments are however in order. Although 
the pairing equation (|40|) does not involve any isotropic 
order parameter, the extended-s wave order parameter 
does have the same point-group symmetry as a pure 
s-wave and if it condenses, will tend to induce charge 
modulation. In a real heavy electron system, the effects 
of Coulomb interaction will renormalise the effective cou- 
pling constant for this channel, eliminating this order pa- 
rameter from consideration. Of the remaining cases, 7^ 
corresponds to a q-CDW order and 7^'^ can be associ- 
ated with spontaneous orbital or line currents between 
the Uranium atoms, as we shall now show. 

Let us consider the current 



let sr^ ( t 



(43) 



from i to j along bond i — j. Orbital order corresponds to 
a non-vanishing circulation of the current in a plaquctte: 

1 



Jc = -(^j.dl^O 



and, therefore, is of the form (Fig. El(^) )• 

^C(X) = i [il2 + hi + j34 + i4l] , 



(44) 



where X) is the position of the centre of the plaquette, 
and the indices (1 — 4) label the corners of the plaquette, 
taking the sense of rotation to be anti-clockwise. 

Now consider the evaluation of the bond variable 
X]o-('^o-(^ + a/2)ccr(x — a/2)). Taking the Fourier trans- 
form each electron field, we obtain 

5](4(x + a/2)c^(x-a/2)) 
= ^ E (cLck'.)e^[''-(^+^/^)-'''-(^-'^/^)] 

k,k',(T 

p 

= ^ E('^(k+k')/2(k - k'))e^[(k-k').x+f^.a] 
k,k' 

= ^^'^"'^ Y.^Pvm^''-^ + (Q - -Q) (45) 

Q) 



where we assumed (pk(q)} = '5k,Q(/0k(Q)) + (Q < 
The current along a given bond is therefore 

j(x-Ha/2,x-a/2) 
= e»'^-^^(pk(Q))2sin(k.a) + (H.c) (46) 

k 

Averaging the currents anti-clockwise around a plaquette 
centred at X, we arrive at 



/c(X) =/cexp[iQ.X] +H.C 



(47) 



where 
^c(X) 



. et 

et 
Nh 



^(/Ok(Q)}[sx sm{Qya/2) - Sy sin{Qxa/2) 

k 

^(pk(Q))K7^ + «-7^], (48) 



where we have used the notation Sx.y = siii{kx,ya) and 
a± = ^ [sin{Qya/2) ± sm{Qxa/2)] 

(notice the ordering of the y and x terms). The form 
factor for orbital current order is thus a weighted mixture 
of 7^ and 7^. Using Eq.llSHl) to simplify Eq. (gHJ, we 
obtain a relation between the orbital current and gap. 



/ = 



eAc i_ 



(49) 



where 



In actual fact, the relative weight of the two channels in 
the orbital antiferromagnet is not an adjustable param- 
eter. If we calculate the divergence of the current at a 
given node in the lattice, we find that 

V-j(r) 
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j{r + ax, r) - jir - ax, r) + j(r + ay, r) - j{r + ay, r) 
|^(a+A3-a_A4) =0 (50) 



so the choice of Q vector determines the mix of 7'^ and 
7^ symmetry in the orbital antiferromagnet. 

In an itinerant model the condition for instability into 
the hidden order phase will be given by the Stoner cri- 
terion already discussed {JiXoip = 1) where xoip ~ 7 so 
that typically at the transition / — where P 0(1) 
is a constant; this is the relation we used for the cur- 
rent in our earlier phenomenological treatment. The 
form factor 7^ = cos{kx) — cos{ky) corresponds to a 
quadrupolar charge density wave (q-CDW). The par- 
ticular details of the conduction electron spectrum ek 
determine which order parameter has a higher critical 
temperature. For instance, if {Qx,Qy) = (t'', tt), and 
Ek = —t(cos{k,j;)+cos{ky)), which corresponds to a nested 
Fermi surface, then from Equation H41|l the relation for 
To is 



Ji 

8ttH 



Tr-\ky\ tanh 



V2T- 



n+\ky\ 



El. 



.(51) 



Here we used i?k(A = 0) = ek+Q/2 — i(sin(fc^)-f sin(fcy)). 
In this particular limit, we can explicitly verify that or- 
bital antiferromagnetism has a higher Tq than q-CDW. 
In the real material, the spectrum may differ greatly from 
this simple form, which may result in a preference of q- 
CDW over OAFM. 

Our discussion in this section is based on a weak- 
coupling treatment of orbital antiferromagnetism, which 
is technically only valid in the vicinity of a nesting insta- 
bility. Real heavy electron systems involve interactions 
of a size comparable with the band-width, in which the 
vicinity to nesting will no longer be a requirement. Prac- 
tical modelling of these situations will require alternate 
strong-coupling methods, such as methods based on a 
Kondo lattice model. It is however interesting to note 
that that both symmetry and microscopic toy treatments 
appear to point to quadrupolar charge density wave and 
orbital antiferromagnetism as the leading contenders for 
hidden order in URu2Si2- 



V. FLUCTUATIONS AND NESTING 

The sharpness of the phase transition in URu2Si2 indi- 
cates that fluctuations do not make a significant contri- 
bution to thermodynamic properties. From the observed 
specific heat anomaly, the region of fluctuations is cer- 
tainly smaller than AT O.IK, so that tg = AT/To < 
2^. This is an unusual situation in the general con- 
text of local moment magnetism, where broad fluctuation 
regions are generally seen in the speciflc heat anomaly. 
This result is sometimes taken to indicate that the hid- 
den order involves a nested Fermi surface ^^i^ However 
band structure calculations have not revealed any signs 



of a nested Fermi surfacejS and it is difficult to see how 
such a condition might occur naturally in the complex 
band-structure of an f-electron system. Sharp mean-field 
transitions are generally taken as an indication of a large 
coherence length scale associated with fiuctuations. In 
insulating systems (e.g. ferroelectrics) this arises from 
the long-range nature of the interaction. In supercon- 
ductors and in nested charge density wave systems, the 
long coherence length ^0 = i'f/A is a consequence of 
the non-local order parameter response of the itinerant 
electron fiuid. 

So can the sharpness of the specific heat transition can 
be used in URu2Si2 be used to infer the presence of nest- 
ing URu2Si2'^. In fact, as we shall now see, a careful ex- 
amination of the Ginzburg criterion for this system shows 
that while we may confirm that the ordering is itinerant 
in nature, the small size of the heavy electron Fermi en- 
ergy means that we do not need to invoke nesting to 
understand that sharpness of the transition. 

The Ginzburg criterion for a phase transition is given 

by 



to 



1 



[(eo/a)'^((55/fcB)]2/(4-'i) ' 



(52) 



or in three dimensions, 



{£.olaf{5S/kB?' 



[d = 3) (53) 



Here a is the lattice spacing, 5S is the entropy associ- 
ated with the phase transition and ^0 is the coherence 
length of the order parameter. Microscopically, ^0 is 
determined from the Gaussian fluctuation term in the 
order-parameter expansion of the free energy. 



AF 



1 



(54) 



where t = 4f — 1 and a is a normalisation constant. The 

Jo 

Gaussian coefficient in the integral is directly related to 
the static susceptibility of the order parameter 



(55) 



The relationship between the coherence length xo find 
microscopic quantities depends markedly on the under- 
lying physics. In insulators, ^0 tends to be determined 
by the range of interaction of the order parameter, but 
in itinerant systems, it is determined by the non-local 
order-parameter polarisation that develops in the elec- 
tron fluid. 

For example, in a local-moment antiferromagnet, with 
interaction H ^ \ SqS^q, 

Xq'=t^B'iT + J,) (56) 
When we expand around the unstable q vector, q — Qq 



J{q) = 0cil 



(57) 
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where Oc = —To is the Curie constant and kT^ the effec- 
tive range of the interaction. With this form, we see that 
for insulating systems, the coherence length = k,^^ 
becomes the range of the interaction. For short-range in- 
teractions, this reason, the breadth of fluctuation region 
is generally large. In insulating systems, narrow fluc- 
tuation regimes are therefore associated with long-range 
interactions. 

By contrast, in itinerant electron systems the order- 
parameter susceptibility generally takes the form 



(58) 



where g is the strength of short-range interaction be- 
tween electrons in the channel corresponding to the or- 
der parameter and xoi> takes the form given in (|42|l . It 
is the momentum dependence of Xo{<l) that determines 
the Ginzburg criterion in itinerant systems. To under- 
stand the role of nesting, let us consider a Fermi surface 
in which the departure from dispersion is measured by an 
energy scale ^, (e.g ek = —2t{coskx -\- cos ky) — /i) then 
the dispersion satisfies 

-Ek-Q = Ek + (59) 

so that the bare susceptibility 1)42(1 is given by 

xo,(q)^f:i7k%/.p (^o) 

^ ek+q/2 - £k-q/2 + 

For /i = and q = 0, this integral is logarithmically diver- 
gent at T = 0, and given by xotpiQ = 0) ~ In (^) 
at finite temperatures. Finite q modifies the Fermi func- 
tions in 1(60(1 . so that 



so that 



'.9 + Xo^ 



HD/To) 



t + 



(61) 



(62) 



and the free energy expansion takes the form 



) (63) 



so by comparing with ((54(1 . we see that for a nested sys- 
tem, the coherence length takes the "BCS" form 



To' 



(64) 



When >> To, then we must replace To ImI i'^ the 
Landau-Ginzburg expansion, i.e. 



AF 





(65) 



from which we see that the coherence length is given by 



Co 



VF 



nestedOj 



(66) 



where we have replaced jjjj- ~ a, so loosely speaking, the 
absence of nesting replaces the coherence length by the 
geometric mean of the BCS coherence length vf/Tq and 
the lattice spacing. 

Let us now return to our case, URu2Si2- Here, using 
the three dimensional form of the Ginzburg criterion, and 
taking URu2Si2, SS ~ O.I/cb, so that 



to 



100 



(67) 



Suppose the fiuctuation region is less than O.IK, i.e. 
to < {0.1K/20K) - 1/200, then a lower bound for the 
coherence length is 

eo/a~(iG/100)-^-(2xlO'^)i/6~5 

Clearly, the presence of the sixth power in the Ginzburg 
criterion, means that only modest coherence length is 
required to account for experiments. Were the hidden 
order strictly associated with the local moments, then 
we would expect (.o/a ~ 1, and clearly, the absence of 
fiuctuations is sufficient to rule this case out. The high 
pressure magnetic phase transition does in fact show clear 
signs of Ising fiuctuations, and in this region, it would ap- 
pear that the ordering transition is indeed local in nature. 
However, we can account for the coherence length of the 
hidden order transition by appealing to itineracy, without 
nesting. By assuming that u_f/A ~ « 25a. Taking 
eF ^ 10'^ K, consistent with the heavy mass m* /nie ^ 60 
and A ~ lO^i^, we are clearly in the right range. From 
these arguments, we see that a correlation length of order 
5 lattice spacings is fully consistent with a system that is 
un-nested but itinerant. We conclude that the sharpness 
of the hidden order phase transition in URu2Si2 only im- 
plies itineracy. Indeed, there are a number of heavy elec- 
tron systems with sharp thermodynamic transitions and 
commensurate magnetic order, indicative of un-nested 
Fermi surfaces, such as 1/2211,11^ and UPd2Al3& In 
each of these cases, it is most likely the itineracy alone 
that is responsible for the narrow fiuctuation regime. 



VI. DISCUSSION 

The observed Fermi liquid behaviour for T > Tq, the 
sharp nature of the transition and the large entropy 
loss point to the hidden order as a general density-wave 
with itinerant excitations formed from the local spin and 
orbital degrees of freedom of the uranium ions and f- 
electrons. Motivated by nuclear magnetic resonance mea- 
surements, we have expanded on our proposal (with J. A. 
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Mydosh) of the hidden order as incommensurate orbital 
antiferromagnetism and have provided technical details 
for our predictions for elastic neutron scattering. Next 
we have turned to a microscopic description of the hidden 
order. After discussing symmetries and allowed particle- 
hole pairings in general terms, we studied the develop- 
ing of these ordering in the setting of a toy single-band 
t— J model within a weak-coupling approach. Within this 
framework, selection between q-CDW and OAFM order- 
ing is not possible, though the situation may be different 
in the (experimentally relevant) strong-coupling regime. 
As discussed in Sec lIII Bl density wave instabilities such 
as q-CDW and OAFM can account for the large entropy 
loss observed at the transition {SS w k%ToN*{0)) if the 
density of states at the Fermi surface, A^*(0), is large 
(as is the case in a heavy Fermi-liquid), and there is a 
substantial gapping of the Fermi surface. 

The weak-coupling model we considered requires the 
nesting of a significant part of the Fermi surface. This 
requirement can be relaxed if the coupling is strong. In- 
deed it seems that a strong coupling description might be 
more appropriate for URu2Si2, since the transition tem- 
perature To is an order of magnitude smaller than the gap 
A unlike a weak coupling description where Tq is more 
comparable with A. Unfortunately here it is difficult to 
perform controlled calculations in this regime, and thus 
experiment is crucial for discerning between these two 
competing scenarios of quadrupolar charge density wave 
order and orbital antiferromagnetism. In Sec. II we stud- 
ied the consequences of OAFM for the neutron scattering 
structure- factor S(q) and NMR at the Si and Ru sites. 
No particular microscopic model was assumed here, so 
the analysis is applicable for any coupling. NMR obser- 
vations were used with our OAFM model to predict an 
incommensurate wavevector for orbital ordering which 
may be verified by neutron scattering measurements. To 
date, these predictions remain untested, as current ex- 
perimental resolution is insufhcient to observe the antic- 
ipated signal level from an OAFM42i^ Here we identify 
a region of momentum space where elastic neutron scat- 
tering probes will clearly be able to distinguish between 
a spin density wave and an OAFM with current signal- 
noise levels. This prediction for orbital anferromagnetism 
remains a challenge for future experiments. 

Our proposal of orbital antiferromagnetism is strongly 
motivated by the inhomogeneous line-broadening ob- 
served in ambient pressure NMR^^ and there are ques- 
tions associated with this experiment that concern us 
greatly. In particular, the local fields measured via NMR 
in epoxied powdered samples are an order of magnitude 
larger than those probed by muon spin resonance or nu- 
clear magnetic resonance in single-crystal ones. One in- 
teresting possibility is "motional narrowing" . The pro- 
posed orbital antiferromagnetic order is incommensurate 
and quite similar in its current patterns to a flux lattice 
of core-less vortices, where the absence of vortex cores 
weakens the pinning effect of disorder. In single do- 
main crystals an incommensurate orbital antiferromag- 



net should then be weakly pinned, giving rise to large 
thermal motion.-^ The probed local-fields will then be 
"motionally narrowed" , i.e their time-average will be sig- 
nificantly reduced in magnitude relative to their its static 
counterpart. One of the predictions of this scenario, is 
that the the muon or NMR linewidth will increase sys- 
tematically with disorder - an effect that might be tested 
using radiation damage to systematically tune the disor- 
der in a single sample. 

Since we started working on this project, there have 
been a number of new experiments which may place 
further constraints on the nature of the hidden order 
in URu2Si2- In particular, recent magnetotransport 
measurements^^ indicate an unusually large Nernst sig- 
nal in URu2Si2 that develops at T = Tq. This kind of 
behaviour has also been seen in the pseudogap phase of 
underdoped cuprate superconductors. In the case of the 
cuprate superconductors, this is most likely an effect of 
the Magnus force on the pre-formed pairs in the pseudo- 
gap. However, here, the absence of any superconductivity 
makes it far more likely that the giant Nernst effect seen 
here is a property of the quasiparticles in the presence of 
the hidden order parameter. These new results clearly 
place an important constraint on the microscopic nature 
of the order parameter. 

Recent high magnetic field studies^ have raised ad- 
ditional questions about the hidden order in URu2Si2- 
Application of high magnetic fields confirms that the hid- 
den order persists to significantly higher values than does 
the remnant antiferromagnetism, affirming the two-phase 
scenarioiiS Moreover the application of still higher fields 
leads to a profusion of new hidden order phases that may 
well cloak a field-induced quantum critical point. At the 
current time, it is not yet clear whether the proposed 
quantum critical point is a consequence of the loss of 
hidden order, or whether it might arise from the close 
vicinity to a quantum critical end point (as is the caseS^ 
with SrRu204.) 

The hidden order mystery in Uranium Ruthcnium-2 
Silicon- 2 can be regarded as part of a much broader set 
of longstanding problems that our community faces in 
the context of highly correlated materials. Coexisting 
forms of hidden order, novel metallic states manifested 
by unusual resistance and magnetotransport properties, 
field-induced quantum phase transitions- each of these 
features, present in URu2Si2, are manifest themselves 
in a wide range of other strongly correlated materials, 
such as the cuprate superconductors, strontium ruthen- 
ate, magneto-resistance materials, and many other heavy 
electron systems. URu2Si2 offers an alternative perspec- 
tive on these problems, and optimistically, its ultimate 
solution will provide part of the key to understanding 
these broader questions. 

Appendix A: SPATIAL DISTRIBUTION OF 
VECTOR POTENTIAL DUE TO OAFM 
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Here we calculate the vector potential A(x) due to orbital order by summing up contributions from currents in all 
links. Consider first the contribution A^^ defined in Eq.JSJ for links along (12). Performing the integral over w gives 



sinh ^ 



a/2 + Xj -X 



sinh 



^^{y-Y,+a/2)^ + {z-Z,)\ 
-a/2 + X.j-x 



^{y-Y,+a/2Y + {z-Z,f 



(Al) 



where we have used the notation Xj = {Xj, Yj, Zj) to denote the co-ordinates of the centre of the plaquette. For the 
link (43) shown in Fig0|we have 



r ^ — ^ 



sinh 



3 

-1 



sinh 



a/2 -^Xj-x 



^^{y-Y,~a/2f + {z^Z,)\ 
-a/2 + Xj - X 



^iy-Y,-a/2r + (z-Z,y 



(A2) 



The X component of the vector potential is then Ax{x) = yl^^(x) + A^'^(x). Similarly the vector potential in the links 
(14) and (23), 



A^^x) 



sinh ^ 



sinh ^ 



a/2 + Y,-y 



^{x-Xj+a/2)^ + {z-Zj)^^ 
a/2 + Yj-y 



r ^ — ^ 



y(.T-Xj+a/2)2 + (z-Z,)2^ 

a/2 + Yj~y 



(A3) 



sinh 



sinh 



^^{x-Xj-a/2Y + {z- ZjY ^ 
-a/2 + X-jij - y 



^{x-Xj-a/2y + {z-Zjy 



(A4) 



yield the y component of the vector potential Ay (x) = A^"^ (x) + A'^'^ (x) . The magnetic field follows straightforwardly 

from B = V X A. 
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